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ABSTRACT 

We investigate a model for the high-ionization state gas associated with observed 
damped Lyman-ct systems, based on a semi-analytic model of galaxy formation set 
within the paradigm of hierarchical structure formation. In our model, the hot gas in 
halos and sub-halos is assumed to be in a multi-phase medium which gives rise to CIV 
absorption, while the low-ionization state gas is associated with the cold gas in galaxies. 
The model matches the distribution of CIV column densities if we assume that the hot 
gas has a mean metallicity log C/H =—1.5, which is the observed mean metallicity of 
damped systems. The same model then leads naturally to kinematic properties that 
are in good agreement with the data, for both the low- and high-ionization state gas. 

We examine the contribution of both hot and cold gas to sub-damped systems 
{Nhi > 4 X 10^®cm“^) and suggest that the properties of these systems can be used 
as an important test of the model. We expect that sub-DLA systems will generally 
be composed of a single gas disk and thus predict that they should have markedly 
different kinematics than the damped systems. We also find that the frequency of 
absorbers drops dramatically for column densities below 4 x 10^® cm“^. These results 
are a consequence of our model for damped Lyman-a systems and we believe they are 
a generic prediction of multi-component models. 

Finally, we find that hot halo gas produces less than one third of Lyman limit 
systems at redshift three. We model the contribution of mini-halos (halos with virial 
velocities < 35kms“^) to Lyman limit systems and find that they may contain as 
much gas as is observed in these systems. However, if we adopt realistic models of 
the gas density distribution we find that these systems are not a significant source of 
Lyman limit absorption. Instead we suggest that uncollapsed gas outside of virialized 
halos is responsible for most of the Lyman limit systems at high redshift. 
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1 INTRODUCTION 

Semi-analytic models of galaxy formation, set within the hi¬ 
erarchical paradigm of structure formation generic to Cold 
Dark Matter (CDM) type models, have been very successful 
in accounting for many of the observed properties of galax¬ 
ies. Previous investigations have focused on the emission 
properties of galaxies in optical, infrared, sub-mm and radio 
wavebands. Some of these successes include matching the 
Tully-Fisher relationship, the optical luminosity function, 
the HI mass function, the galaxy correlation function, and 


distant galaxies (most recentlv Somerville & Primack 1999i 

Cole et al. 

bOOd; iDevriendt & Guiderdoni tOOOl; iNagashima 

et al. 2001 

; Somerville et al. 2001; Springel et alj 2001). 


H owever, with the notable exception of Kauffmann 
( 1996|), these models have not been used to address obser¬ 


vations of absorption systems in much detail, though sev¬ 
eral papers have performed the basic consistency check that 
the total amount of cold gas in these models is at least as 
much as that observed in damped Lyman -g (DLA) systems 
(Baugh et al. 1998; Somerville et al. [2001 ). The importance 


of utilizing absorption system data to confront galaxy for¬ 
mation models should not be underestimated. The optical 
properties of galaxies are determined by a sequence of poorly 
understood physics: star formation, the initial mass func¬ 
tion, supernova feedback, metal enrichment, dust, etc. In 
contrast, absorption systems present an orthogonal means 
of observing the high redshift Universe, and are subject to 
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different theoretical uncertainties as well as different obser¬ 
vational selection effects. Demanding that a single theory 
should be able to account simultaneously for the observed 
properties of 6oth absorption and emission systems therefore 
presents a stringent test for any theory of galaxy formation. 

In this paper, we use a model that has previously been 
shown to successfully account for the pr operties of nearby 
and 2 3 g alaxies identified in em ission (Somerville & Pri¬ 

mack 1 )99; Somerville et al, 2001) to develop a parallel pic¬ 


ture that can also account for the observed properties of 
absorption systems. This paper is the second in a series 


based on this approach. In the first paper (Mailer et al 


basic ingredients of our model. In Section 4 we compare 
our model with the observations. In Section p, we investigate 
the properties and origin of LL and sub-DLA systems in our 
model and identify future observational tests of the picture 
we have proposed. We close with our conclusions in Section 

§ 


2 BACKGROUND 

Early on it was established that damped absorption sy stems 
are an important probe of galaxies at high redshift (Wolfe 


2001, hereafter paper I), we demonstrated that the kine- 


et al, 1986). The numerous subsequent studies of DLAS 


matic properties of DLA systems (absorption systems with have focussed on their frequency, metallicity, and kinemat- 
HI column densities > 2 x 10^° cm~^) as traced by low- 


) Prochaska & Wolfe 

1997, 

1998) could by 

Peroux et al 

(pool), and 

Prochaska Sz Wolf^ (2002 


be understood in this context if the gaseous disks of galax¬ 
ies are very extended, so that a significant fraction of the 
absorption systems arise from lines of sight passing through 
multiple disks. We showed that the metallicity and number 
density of DLA absorbers in this model were also in general 
agreement with observations. In this paper, we turn our at¬ 
tention to the kinematics of the associated high-ionization 
state gas and its correlation with the cold neutral gas. The 
study of the relationship between high- and low-i onization 


state gas in DLA systems has been pioneered by Wolfe & 
Prochaska (2000a q), who found that none of the models 


they explored could successfully match their observations. 

The basic premise of our model is to associate the highly 
ionized gas that gives rise to CIV absorption with the hot gas 
in dark matter halos and sub-halos. Locally, high-ionization 
state gas like CIV is k nown to be associated with galactic 


halos (Chen et al. 2001). Unfortunately, modelling such sys¬ 


tems is rather complicated b ecause the gas is likely to b e 


multi-phase (see for example Mo & Miralda-Escude 1996). 


We do not attempt to model the multi-phase medium in any 
detail, but rather to ascertain if the amount of hot gas in the 
galaxy formation model, distributed in a reasonable manner, 
can produce the observed high-ion kinematics. We find that 
such a model is fairly successful, suggesting that the general 
picture is compatible with the combined constraints from 
the high-ion and low-ion data. 

We then explore the contribution of hot and cold gas 
in our model to lower column density HI systems. We find 
that cold neutral gas produces sub-DLA systems (column 
densities of ~ 4 x 10^® — 2 x 10^° cm“^) and we demonstrate 
how the kinematics of these systems may differ from the 
higher column density DLAS. We also find in our model 
that CIV systems with associated low ion absorption have 
markedly different kinematic properties than those which do 
not, i.e., kinematic information discriminates between DLA 
and sub-DLA systems. 

Lastly, we find that hot halo gas can produce no more 
than one third of observed Lyman limit (LL) systems, im¬ 
plying that the majority of these system s are produced b y 
either uncollapsed gas or mini-halos (e.g. Abel & Mo 1998). 
We model the contribution of mini-halos and find that when 
we adopt a realistic gas distribution and include an ionizing 
background, we find that these halos are unlikely to be the 
source of a large fraction of LL absorbers. 

The paper is organized as follows. Section ^ reviews 
the observations of DLA systems. Section W describes the 


presented the most recent study of their metallicity. Most 
of the kinematical data has been obtained by Wolfe and 
Prochask a in an amb itious program using the HIRES spec¬ 
trograph ( Vogt 1992 ) on the Keck Telescope. They acquired 
velocity profiles for a number of metal lines associated with 
each DLA system. They found that low ionization state lines 
(e.g. Sill, AlII, Fell) traced one another well and therefor e 
presumably also the cold HI gas (Prochaska & Wolfe 1996). 
Higher ionization state lines (e.g. SilV, CIV) showed dif¬ 
ferent but correlated behavior, im plying that the y traced a 
different phase of the ga s (see also 


Lu et al 


1996). 


Troctiaska & Wolte (1997) introduced four statistics to 
characterize the gas kinematics: Av, the width containing 
90% of the optical depth; /mm, the distance between the 
mean and the median of the optical depth profile; fedg, 
the distance between the highest peak and the mean, and 
/zpfe, the distance of the second highest peak to the mean 
of the profile. The last three statistics are suitably normal¬ 
ized to haye ranges between 0 and 1 or —1 and 1. They also 
introduced three statistics to quantify the correlations be¬ 
tween the high- and low-ion kinematics (Wolfe & Prochaska 
2000a). These are 


• 5v, the offset between the mean velocity of the high- 
and low-ionization state gas, 

• /ratio, the ratio of the velocity widths of the high- and 
low-ionization state gas, and 

• fhe cross-correlation between the two gas states. 

Previou s semi-analy t ic mo delling within the CDM cos¬ 
mogony by Kauffmann (1996) suggested that the major¬ 
ity of the cross section to DLA systems should come from 
small m ass halos with virial veloc ities between 35 and 50 
kms”'^. Prochaska & Wolfe (1997) investigated a number 
of models to explain the kinematics of the low ionization 
state gas. They concluded that th e distribution of veloci¬ 
ties predicted by Kauffmann (1996) was strongly ruled out 
by the data. The only model they investigated that was 
compatible with the data was one in which the absorp¬ 
tion systems arose f rom a thickened, rapidl y rota ting disk 
(Kot — 200kms“^). ledamzik fc Prochaska] ( 199^ ) demon¬ 
strated that most CDM-based models for galactic disks were 
incompatible with the observed velocity widths of the DLA 
systems — under the assumption that a single disk gave 
rise to each absorption system. Other workers, however, 
found that multiple component models, in which some of 
the DLAS were produced by lines of sight passing through 
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multiple disks or proto-galactic clumps, could explain the 
low-ion observations within the context of hierarchical struc¬ 
ture fo rmation (pape r I; 


Mirald^ ^Escud^ 


199' 


Haehnelt et al 


199f; 


McDonald & 


Gardner et ah' (pOOl) used cosmological hydrodynami- 


cal simulations to determine the cross section of DLA and LL 
systems. They found their simulations could produce the ob¬ 
served number density of DLA and LL systems in a ACDM 
cosmology if sufficiently small mass halos were considered 
(60kms“^ for DLA systems and 30kms“^ for LL systems). 
They found the cross section as a function of halo virial ve¬ 
locity co uld be fit by a powe r law with a slope of « 1.6. 
However, Prochaska & Wolfe (2001) have pointed out that 
combining this distribution of cross-sections with the veloc¬ 
ity w i dth to virial velocity relationship found by Haehnelt 
et al. ( 1998 [) and in paper 1 is not compatible with the kine- 


distribution of velocity widths (Haehnelt et al. 2000) which 
is consistent with what we f ound for our m odel in paper I. 
The different value found by 


Gardner et al 


is probably due 
to the fact that there is effectively no feedback in their sim¬ 
ulations leading to much higher baryon fractions in small 
mass halos. In models with efficient feedback we would ex¬ 
pect the contribution from these low mass halos to be sub- 
stantially decreased. _ 

Woite & Froctiaska! ( 2000b| ) found that none of their 
models including the thick disk model were consistent with 
the high-ionization state gas kinematics. They suggested 
that the multiple component model might be more success¬ 
ful, which we now explore. 

Another long-standing question pertaining to QSO ab¬ 
sorption systems is their connection with the galaxy popula¬ 
tion identified in emission. At 2^1, there is compelling ob¬ 
servational evidence th at metal line and LL absorbers are as - 
sociated with galaxies ( [Steidel et al. 1997; Chen et al, 2001). 
A theor e tical model has been proposed by |Mo fc Mirald^ 
Escude ( 1996 ) wherein these systems are produced by gas 
clouds in the halos of galaxies. At higher redshifts, it has 
been propo sed that many o f the absorbers may reside in 
mini-halos (Abel & Mo 1998). We also investigate the origin 
of lower column density systems in the context of our model. 


is labelled as either “hot” (at the halo virial temperature) 
or “cold” (T^IO** K). 

When halos merge, the central galaxy of the largest pro¬ 
genitor halo becomes the new central galaxy, and all other 
galaxies become satellites. Satellites lose angular momentum 
due to dynamical friction and gradually fall towards the cen¬ 
ter of the halo, where they may eventually merge with the 
central object. 

Cold gas is converted to stars using a simple empirical 
recipe, at a rate proportional to the total mass of cold gas 
in the galaxy. Several d iffere nt rec ipes were investigated i n 
Somervifle & Primack (1999) and Somervifle et al, (2001), 


but here we adopt the “collisional starburst” recipe, which 
was shown to produce the best agreement with the observed 
galaxy population at z ~ 3, and with the total mass of neu- 


matic c ata. A slope of 2.5 is needed to produce the correct _ et al. 


tral hydrog e n implied by observat ions of DLAS (Somerville 


2001 


Primack et al 


2001). In this model, cold gas 


in isolated disks is assumed to form stars with a low effi¬ 
ciency (simifar to that in nearby observed spiral galaxies), 
while a more efficient starburst mode is triggered by galaxy 
mergers. Gas that has settled into a disk may be heated and 
removed by supernova (SN) feedback. The efficiency of the 
SN feedback is assumed to be inversely proportional to the 
potential well depth of the galaxy. Each generation of star 
formation produces a certain yield of heavy elements, which 
are mixed with the cold gas or ejected by SN into the hot 
gas halo. 

After reionization, cooling only takes place in dark mat¬ 
ter halos with virial velocities > 40 kms“^|^ because of 
suppression of gas collapse and cooling by the UV back¬ 
ground. 

All of our analysis has been done assuming the currently 
favored ACDM cosmology (e.g. Primack 2002) with Ha = 
0.7, Hm = 0.3, Hb = 0.038,/i = 0.7, and as = 1.0 and at 
a redshift of three. The values of the free parameters that 
control the efficiency of star formation and SN feedback. 


the chemical yield, etc., are as adopted in Somerville et al 


(^001 ) based on low redshift galaxy observations. 


Below we describe in some detail how we incorporate a 
model for both the low- and high-ionization state gas in ab¬ 
sorption systems into this picture. A schematic illustration 
of the model is shown in Figure 1. 


3 THE GALAXY FORMATION MODEL 


For our analysis we use the se mi-an alytic model desc ri bed i n 
Somerville & Primack (1999) and ^omerville et al, (2001). 
The backbone of the model is a Monte Carlo realization of 
a “merger tree”, which represents the buil d-up of halos over 
time through merging and mass accretion Somerville & Ko- 


latf (1999). Initially, the hot gas is assumed to be distributed 
like the dark matter (here, a singular isothermal sphere), to 
be uniformly at the virial temperature of the halo, and not 
to have any substructure. The cooling radius (the radius 
within which gas is dense enough to have cooled) is then 
calculated using the radiative cooling function for atomic 
gas in collisional equilibrium. When a halo contains more 
than one galaxy, gas is ass umed to cool only onto the cen¬ 
tral o bject. Recent studies ( Toshida et al. ^002| ; [Helly et al 


2002) suggest that this simple model is in good agreement 
with cosmological hydrodynamical simulations, at least in 
the absence of feedback. All gas in these simplihed models 


3.1 Low-Ionization State Gas 

We assume that the low-ionization state absorption is asso¬ 
ciated with the “cold” gas component of the galaxy forma¬ 
tion model. In paper I, we explored a number of models for 
the distribution of cold gas in each galaxy and found that 
only by having very extensive gas disks was it possible to 
produce enough overlapping cross-section for a line of sight 
to a distant quasar to pass through multiple disks and thus 
produce low-ion kinematics in agreement with the observa¬ 
tions. We found that a model where the cold gas was dis¬ 
tributed such that its surface density decreases proportion¬ 
ally to R~^ until a truncation radius Rt of column density 


^ The virial velocity is just the circular velocity at the virial ra¬ 
dius, which for isothermal spheres is the circular velocity every¬ 
where. 
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Figure 1. Gas in a typical halo at 2 = 3 seen in projection. The 
filled ellipses are cold neutral gas disks, while the solid circles 
mark the edge of the spherically distributed hot gas in subhalos. 
The entire halo is also filled with hot gas out to the virial radius 
(the dashed line). Three lines of sight through the halo illustrate 
how different components can give rise to DLA, sub-DLA and LL 
systems in our model (see section m. 


A^t = 4 X 10^® cm“^, produced absorption sy stems consis- 
tent with th e low ion kinematics observed by Prochaska & 
Wolfe (1997, 1998) as well as with the observed f{N) dis- 
tribution^ and metallicities of DLA systems (see paper I for 
more details). 

We take the vertical (z) distribution of the cold gas to 
be exponential with a scale length Zs such that the disks are 
thin. Thus the gas distribution for each galaxy is given by 


{R<Rt). 


( 1 ) 


We distribute the cold gas in 10 clouds to properly sample 
the velocity held. Each cloud is given a small, ID random 
velocity of aiow = 10kms“^. 

We refer to the cold gas as being in “disks”, and we im¬ 
plicitly assume that the cold phase is concentrated and ro- 
tationally supported when implementing the star formation 
recipes. However, for the purposes of the kinematic mod¬ 
elling, the cold gas need not be rotationally supported — 
by disks we simply mean battened systems. We emphasize 
that this description is a crude simplification of the distri¬ 
bution and dynamics of cold gas in real galaxies, especially 
in interacting galaxies, which show rather complex behavior 


^ f{R) the number of absorption systems with column densities 
between N and N+dN found per comoving path length. 


(Hibbard 2000). In the spirit of semi-analytic models, it is 
our hope that such a simplified description can still capture 
the general behavior of the cold gas when averaged over a 
large number of systems. 


3.2 High-Ionization State Gas 


The high-ionization state absorption is assumed to arise 
from the hot gas in halos. In the standard galaxy formation 
model, all the hot gas is assumed to be at a single tem¬ 
perature and only a single gas phase is considered. These 
simplifications may be reasonable for calculating the cool¬ 
ing rate of the gas, but for describing kinematics we will 
need to consider the substructure of the hot gas, and we 
will assume that the majority of it is in a phase suitable 
for CIV absorption, in general at temperatures less than the 
virial temperature of the halo. We will continue to refer to 
this as hot gas because it is dynamically hot even if parts of 
it have cooled. 

To compare to the observations we must convert the 
hot gas mass to a CIV column density, Nciv- The column 
density of CIV along a line of sight is related to the column 
density of hydrogen by 

Nciv = fcivZhgNn, (2) 


where Zhg is the metallicity of the gas and fciv is the frac¬ 
tion of the gas in a state that produces CIV absorption. 
The column density in HI is related to the total amount of 
hydrogen by Nhi = (1 — x)Nh, where x is the ionization 
fraction of the gas. The quantities x, fciv and Zhg can vary 
locally in the gas, but we will assume that suitably averaged 
quantities can be defined to give the global relationships we 
use here. 

We emphasize that only the metal lines in the systems 
are measured, which creates an inherit degeneracy between 
the total amount of gas, its metallicity and its ionization 
state. For simplicity we will treat the ionization state and 
metallicity to be uniform within each halo and therefore 
Nciv is directly proportional to Nh- The reader should bear 
in mind however that an increase in Nh can equally well be 
thought of as an increase in fciv or Zhg- We well return to 
this issue in § ^ 

For simplicity we take fciv = 1, noting that a smaller 
value of fciv can be directly offset by an increase in the 
hot gas metallicity. We fix Zhg to be —1.5 dex for all hot 
gas in all halos, which produces a good match to the CIV 
column density dist ri bution in DLA systems as measured in 
Wolfe & Prochaska (2000a). The mean metallicity of DLA 
syste ms at2>2is«—1.5to—1 dex (Prochaska & Wolfe 
2002), which is consistent with our assumed Zhg- However, 
because of the degeneracy mentioned above, this value is 
directly proportional to our choice of fciv- Thus this can 
be thought of as a minimum metallicity since fciv < 1- 
A realistic estimate of fciv might be « 0.35 which would 
then require Zhg ~ —1.0 dex. Alternatively, the hot gas 
metallicity could be lower if there is more hot gas in halos 
than calculated by our simple recipes. 

The galaxy formation code calculates the metallicity of 
the hot gas as part of the chemical enrichment model de¬ 
scribed briefly above. We find that if we use the metallicities 
calculated in this way (which typically show a large scatter 
at fixed halo mass), we obtain too large a spread in the CIV 
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column density distribution. This is not too disturbing, as 
the modelling of the reheating and ejection of metals by 
supernova feedback is extremely uncertain. The galaxy for¬ 
mation model was calibrated to produce the correct mean 
observed stellar, cold gas, and hot gas metallicities in 2 = 0 
galaxies and clusters, but this is no guarantee that it pro¬ 
duces the correct results for these quantities at 2 = 3, or the 
correct dispersion in the mass-metallicity relation. 

However, it is also possible that correlations between 
fciv, Zhg and the gas density conspire to reduce the spread 
in observed CIV column densities. An example of this kind of 
conspiracy would be a starburst which produces high metal- 
licity in the hot gas but also reduces the hot gas column 
density or changes the ionization state of the gas such that 
CIV is no longer the preferred ionization state of carbon. 
Unfortunately, the large number of uncertainties mean that 
few constraints can be placed on any aspect of the mod¬ 
elling that goes into determining CIV column densities. We 
discuss constraints that can be placed on these quantities in 
section 

We will assume that the hot gas is distributed like the 
dark matter, i.e. density proportional to r~^. For our kine¬ 
matic study we would like to consider hot gas associated 
with each galaxy in the halo, while in the galaxy formation 
model, all hot gas is assumed to be associated with the cen¬ 
tral galaxy. To do this we assume that the hot gas in each 
subhalo is related by a parameter fsub to the total amount 
of dark matter in the subhalo. Thus the mass of hot gas, 
rrihg in a given subhalo is 


'b^hg — fsubfT^dm 


Mhg 

Mdm 


(3) 


where rudm is the dark matter mass of that subhalo and 
Mhg and Mdm are the total mass in hot gas and dark matter 
respectively, summing over all subhalos and the parent halo. 
The remaining hot gas is then associated with the parent 
halo. In our realizations we find that subhalos can contain 
up to 2/3 of the halos mass, requiring fsub to have values 
between 0 and 1.5. In the case that fsub ~ 0, there is no hot 
gas associated with the subha los and our high-ion model is 
similar to the one explored in Wolfe & Prochaska (20001). 

We note that this is far from a model of the hydrody¬ 
namics of the hot gas in the halo, which is rather difficult to 
model especially including the poorly understood physics of 
supernova feedback. A successful model of gas in the halo of 
the Milky Way, where there is plentiful data, still eludes re¬ 
searchers (Collins et al 2002), so we must accept that at the 
present moment we can only hope to sketch a rudimentary 
framework of the situation at high redshift. 

For the purpose of simulating spectra we assume that 
a fraction fciv of the hot gas is no longer at the virial 
temperature but has cooled into pressure confined clouds 
and is photoionized with a temperature « 5 x 10“* degK. 
We use 15 clouds to properly sample the velocity field and 
assume clouds have random velocities, aciv, proportional 
to the halo’s circular velocity (see Table ^). 

We compare our model only to the CIV data which 
best samples the hottest phase of the gas. The fact that 
there are s ome differences in the hi gh-ion kinematics be¬ 
tween ions (Prochaska & Wolfe 2002, e.g. SiIV vs. C IV) 
implies variations in the ionization state or metallicity of 


the gas. Future modelling of this may help us to understand 
metallicity, density and ionization gradients in the hot gas. 


3.3 Mini-halo Model 


Mini-halos (Vvir < 35kms“^) are usually not considered 
in discussions of galaxy formation because their exp ected 
luminosities make them difficult to observe. However, Abel 


& Mo (1998) have shown that there may be sufficient cold 
gas in mini-halos at high redshift to comprise the LL systems 
at those redshifts. 

Recently the issue of mini-halos has also been widely ad¬ 
dressed because of what has been named the dwarf satellite 
problem', the fact that N-body simulations expect hundreds 


of m i ni-subhalos in a M ilky Way sized halo (Klypin et al 


1991 ; Moore et al, 1999 ), but only a couple dozen or fewer are 


observed. While many exotic solutions to this problem have 
been proposed, the simplest is that the extragalactic UV 
background suppresses the accretion and cooling of gas in 
these mini-halos, and therefore stars form in only the small 
number that managed to collap se before the epic of reion¬ 


ization. Models of this scena rio (Bullock et al. 


et al. 


2002 


2000 


Benson 


Somerville 2002) have been fairly successful in 


explaining why only a small fraction of mini-subhalos would 
have formed stars. 

To explore the contribution of mini-hal o s to L L systems 
we use the squelching model of ^omervilie (2002) to deter¬ 
mine the amount of hot and cold gas in mini and low-mass 
(35kms“^ > Vvir > 50 kms~^) halos. This model uses a 


fitting function found by Gnedin (2000) in hydrodynamical 
simulations to determine the amount of gas that can accrete 
onto a dark matter halo in the presence of an ionizing field. 

We assume that gas will only contract to density pro¬ 
files steeper than the dissipationless dark matter. Therefore 
to explore the maximal contribution of these halos to LL 
systems, we take their projected density profiles to go as 
R~^ like the cold gas disks and hot gas in galaxies described 
above (equation |^. We truncate the gas at the virial radius 
or we consider cases where the gas is truncated at a specified 
column density, Nt- In the latter case the truncation radius 
can then be evaluated by 


Rt = 


Mg. 


2'KNtmH 

where mn is the mass of the hydrogen atom. 


(4) 


4 RESULTS 

We pass multiple random lines of sight through each halo to 
produce Monte-Carlo realizations of spectra in the manner 
described in paper 1. Enough realizations are performed to 
give us at least 10,000 sight lines with DLA systems in them. 

Figure]^ shows the distribution of CIV column densities 
associated with DLA systems produced in our model. We see 
that the model distribution is in good agreement with the 
observations. We chose the hot gas metallicity (—1.5 dex) so 
that the bottom range of log column densities is about 12.5 
as seen in the data. The highest bin contains all systems 
with log column densities greater than 14.8 as this is where 
the CIV doublet becomes saturated and only lower limits on 
the column density can be determined. Of course we expect 
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log Nciv 


Figure 2. The distribution of CIV column densities for the data 
and the model are shown. The model is normalized to have the 
same number of systems as the data so that Possion errors can 
easily be estimated. The CIV doublet is saturated for systems 
with log column densities higher than 14.8; these systems are 
binned together. We have chosen a hot gas metallicity of —1.5 
dex so that the low-end cutoff of the model matches the data. 
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that there is some spread in the hot gas metallicities. If we 
take Zhg to be distributed as a Gaussian in log metallicity 
then a mean of —1.5 dex and a standard deviation of 0.5 is 
in equally good agreement with the observations. 

We can place a lower limit on the amount of hot gas 
in galactic halos needed to match the Nciv distribution by 
taking the highest reasonable values for fciv and Zhg- Since 
we have already started with the maximum value of fciv, 
the only way to lower the amount of hot gas is to increase 
the metallicity, Zhg- If we believe that Zhg must be less 
than a third solar then we could get the same results with 
one tenth the mass of hot gas. Any model with less than this 
amount of hot gas would have difficulty producing enough 
CIV absorption. 

Figure ^ shows the results for the most important kine¬ 
matic characteristics of the high-ion gas: Aviow, Avciv, Sv, 
and fratio (see Section]^ for definitions of these quantities). 
We show the results for five variants of our model (a-d, and 
d-h), summarized in Table In the first four models we vary 
the parameters aciv and f^ub as detailed in Table All 
four models produce acceptable fits to the kinematic data 
and none of them can be rejected with high confidence. The 
statistical tests are sensitive to complimentary aspects of the 
DLA kinematics and place relatively tight constraints on the 
model parameters. For example, increasing aciv (model b) 
improves the fratio but lowers the value of the Avciv statis¬ 
tic. Making the gas more clumped in subhalos (model d) im¬ 
proves the agreement with the observations but we cannot 
rule out a model with no cl umping (model c). 

Previous attempts by Wolfe & Prochaska (2000b) to 
model the high-ion gas of the DLA failed for scenarios in¬ 
volving a cold, rotating disk surrounded by a hot gas halo. 


Figure 3. The four most discriminatory tests are shown for high 
and low ions. These are the Avi^n,, Avciv, and fratio tests. 
Results are shown for four models where the parameters aciv 
and fsub have been varied (see Table ^1. Also included is a model 
(d-|-) where we have added an additional source of CIV absorption 
associated with the cold gas disk. All the models are consistent 
with the observations. Increasing fsub and the addition of some 
CIV absorption in the gas disk improves the models’ agreement, 
but not drastically. Adjusting other parameters like aciv can 
improve one statistic {Avciv in model b) but at the expense of 
another statistic {/ratio)- 


In short, the rapidly rotating disks favored by the low-ion 
kinematics lead to a bimodal distribution of Sv values which 
is ruled out at a high confidence level by the DLA observa¬ 
tions. The source of this inconsistency is the fact that the 
hot gas velocity field is centered at the rest-frame velocity of 
the halo while the rotating cold gas has a preferred velocity 
of ~ ±Vvirl‘2- This offset between the cold disk gas and hot 
halo gas is not found in the observed Sv distribution. For 
the disks in our model, this offset also applies. The key dif¬ 
ferences between our model and the single disk models with 
respect to this issue are: 

(1) the rotation speeds of the SAM disks are consider¬ 
ably smaller t han the v„ir ~ 200 km/s required by single 
disk scenarios (Prochaska & Wolfe 199^ ); and 

(2) the centroids of the hot and cold gas velocity fields 
converge as the number of disks intersected along a sight 
line increases. 


While a large range of parameters (Table lead to 
acceptable models, the agreement with the statistic is 
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Model 

^CIV 

fsub 


Avciv 

Sy 

fratio 

a 

Hvir / 

1.0 

0.30 

0.75 

0.09 

0.10 

b 

Hyir 

1.0 

0.30 

0.19 

0.08 

0.56 

c 

Hyir 

0.0 

0.31 

0.12 

0.03 

0.24 

d 

Vyij> 

1.5 

0.32 

0.44 

0.07 

0.50 

d+ 

Vyij> 

1.5 

0.30 

0.66 

0.28 

0.31 


Table 1. Five models where the parameters crciv fsub are 
varied. In the d+ model CIV absorption associated with the cold 
gas disk is included. The last four columns give the KS proba¬ 
bility that the model is consistent with the data. All five models 
produce acceptable matches to the data. The d-|- model produces 
significantly better agreement with the statistic than the other 
models. 


marginal. This can be attributed to the relatively large num¬ 
ber of observed DLA with nearly zero offset between the 
mean velocities of the CIV and low-ions compared to the 
models. Figure ^ shows the cross-correlation between high 
and low ions in our models and in the data. We see that the 
four models are in good agreement with the data when the 
velocity difference is greater than 50kms“^ but become less 
correlated than the data for smaller velocities. This implies 
that there are additional correlations between high and low 
ions that we have not including in our modelling. 

We therefore hypothesize that some of the CIV is not 
produced by hot halo gas, but instead is ass ociated with 


the c old gas disk as is seen in the Milky Way (savage et al 


1997). This component would then be highly correlated with 


the low-ion kinematics. We model this component by taking 
its total log column density to be 12.5 and its kinematics to 
be identical to the low ions. We then add this to the CIV 
produced by the hot halo gas. 

Figure ^ shows the cross-correlation for a model with 
this additional component (model d-|-). One sees that the 
agreement with the data for velocity differences less than 
50kms“^ is greatly improved. These models also have more 
success with the Sv statistic (Figure^ than the models with¬ 
out the additional component. Thus we conclude that our 
model without any additions is a good model of the kinemat¬ 
ics for scales larger than 50kms“^. At smaller scales detailed 
modelling of the ISM will be necessary to truly model the 
contribution of high-ion gas correctly. 

Therefore we conclude that the basic picture presented 
here for the origin of the high- and low- ionization state 
gas in absorption systems successfully reproduces all known 
observable properties of these systems (aside from the possi¬ 
ble issue of hot gas metallicities). Next we turn to modelling 
lower column density absorption systems, about which much 
less is known. Here, future observations will be able to test 
the validity of the scenario we have presented for DLA sys¬ 
tems. 


5 SUB-DLA AND LL SYSTEMS 

We have investigated the properties of DLA systems and 
the CIV absorption associated with them. Our model also 
makes predictions relating to gas that is below the DLA col¬ 
umn density limit of log(Aim) = 20.3. It has been observed 
locally that gas disks are truncated at column densities of a 



Figure 4. The cross-correlation, C(r)), between high- and low- 
ionization state gas. The solid line is for the data while the other 
line types are for the five models listed in Table ^ We see that 
the first four models show good agreement with the data down to 
velocities of 50kms“^ but then become less correlated than the 
data. By adding a high-ion component that rotates with the low- 
ion gas (model d+), we increase the correlation at small velocity 
differences, obtaining a better fit to the observed (^{v) and leading 
to an improvement in the (5„ test. 


few xlO^^cm which has been attri buted to photoiniza- 
tion by extra-galactic UV radiation (Corbelli & Salpetei 


1993; Maloney 1993). We have found that when we trun¬ 
cate our disks at 4 x 10'^®cm“^, our model produces good 
agreement with the kinematic data at z ~ 3. One would 
expect the column density at which gaseous disks are trun¬ 
cated to decrease from 2 = 3 to the present as the intensity 
of the UV background drops. 

Systems with log column densities less than 20.3 but 
greater than 19.0 ha ve been referred to as sub-DLA systems 
( Peroux et al. 2001 ). We suggest that sub-DLAS should re¬ 
fer to cold neutral systems like DLA and therefore should 
have column densities greater than a few xlO^® cm~^. The 
cutoff is probably a function of redshift as the strength of 
the UV background and/or density of the gas disks varies. 
Measurements of the ion abundances in these systems should 
be able to tell us where the systems make a transition from 
cold neutral gas to photo-ionized gas. For now we will adopt 
a value of 4 x 10^® cm“® as the lower limit of sub-DLA sys¬ 
tems to be consistent with the column density at which we 
truncate our gas disks in our model. 

At lower column densities, LL systems must arise from 
something other than cold neutral gas. Lines of sight through 
galactic halos that fail to intersect cold neutral gas may give 
rise to LL systems (see Figure |^. It is also possible that LL 
systems are produced by uncollapsed gas not in virialized 
halos, or gas in mini-halos. 


5.1 The contribution from mini-halos 

Figure ^ shows the number per unit redshift dN/dz of mini¬ 
halos as a function of v^ir- We see that the cross section 
for encountering a halo with a virial velocity in the range 
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Figure 5. d^N/dzdV for mini-halos as a function of the halo’s 
virial velocity. The solid line is the cross section of the dark mat¬ 
ter. The dashed and dot-dash lines assume, respectively, that the 
hot or cold gas extend to the halo’s virial radius. The hatched area 
assumes the cold gas is truncated at cm“®. The smaller filled 
area assumes the hot gas is 99% ionized. We see that in realistic 
models of mini-halos they do not contribute much cross section 
to LL absorption. The cold gas in low mass halos may account 
for a substantial fraction of sub-DLA systems if it is truncated at 
a column density lower than that of cold gas in galactic halos. 


Figure 6. The column density distribution f{N) at z = 3 for 
D T.A and T,T, sy s tems in our model. Also shown are the data 


of Peroux et al 


(2001). Note that the three points on the left 


are not actiia.l data, nnints h ut an extrapolation of the dotted 
line (see Peroux et al. [2001 ). The dark line shows the column 
density distribution of cold neutral gas. The dashed and dash-dot 
lines show the column density distribution of hot gas assuming 
95% and 99% ionization respectively. One sees that there is a 
large drop in the number of absorbers with log column densities 
between 18.0 and 19.8 in our model. One also sees that hot halo 
gas only generates about a third of the LL systems at z = 3. 


10 — 50kms“^ at z = 3 is 10 times as large as the observed 
dN/dz for LL systems. 

As we have discussed, halos in this velocity range will 
have difficulty in accreting and cooling gas in the presence of 
a UV background. However, mini-halos which collapse be¬ 
fore reionization may contain cold gas. If we assume that the 
cold gas in these systems extends out to the virial radius of 
the halo we obtain dN/dz = 3.7 (the d ashed line in Figur e 
This is consistent with the results of Abel & Mo (1998). 
It requires that the cold gas remain neutral at column den¬ 
sities far below 10'^®cm“®, which is probably unrealistic. If 
instead the cold gas is truncated at 10^® cm~® then we ob¬ 
tain dN/dz = 0.3, a negligible contribution to LL systems. 
In this case, the mini-halos would account for more than half 
of the systems with Nhi > 10^® cm“®, but it is unclear why 
the truncation value for mini-halos should be less than that 
of galactic halos. 

An examination of the hot gas in mini-halos leads to 
similar results. If the hot gas is taken to be neutral and to 
extend out to the virial radius then it has a rather high rate 
of incidence dN/dz = 4.7. However, once again it is difficult 
to understand why this gas would not be highly ionized. 
In hydrodyn amic simulations of mini-halos in an ionizing 
background, Kepner et al. ( 1999 ) found that at maximum 
the HI column density only exceeds the Lyman limit thresh¬ 
old out to a radius less than one ninth the virial radius. If 
we assume that the hot gas is uniformly 99% ionized then 
dN/dz < 0.06. 

Thus it seems that cold gas in mini-halos may con¬ 
tribute to sub-DLA systems but that in realistic scenar¬ 


ios mini-halos p r oduce relatively few LL systems. Note that 
Gardner et al. (2001) suggested that their hydrodynami- 
cal simulations could account for the number of LL ab¬ 
sorbers if they could be continued down to halos with 
Vvir ~ 30kms“^. The reason these results are different from 
ours using semi-analytic models is that the simulations have 
effectively no feedback; there is much more cooling and early 
star formation than in the semi-analytic models. Also they 
see no effect of the UV photoionizing background on ha¬ 


los with Vvi' 
mini-halos. 


> 50 km s but cannot resolve low mass and 


5.2 The contribution from hot gas in galactic 
halos and uncollapsed gas 

The hot gas in halos which produces the high-ion gas in DLA 
and sub-DLA systems will also generate LL systems in lines 
of sight that do not intersect cold gas. The frequency of 
fines of sight that intersect hot gas in our model is only 0.78 
per unit redshift, well below the observed number density of 
2.0 ± 0.5 for LL systems. Furthermore the hot gas is highly 
ionized so that depending on the ionization fraction x these 
systems may not have a high enough column density to be 
LL systems. 


Frochaska (1999) has observed that x = 0.97 ± 0.02 


in a LL system. If we assume that x is within this range 
for the hot gas in halos then the resulting column density 
distribution is shown in Figure |^. The dashed line and the 
dot-dashed line show the distribution taking x of 0.95 and 
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0.99 respectively. Column densit i es hav e only been measured 
for DLA systems; Peroux et al. (2001) has extrapolated the 


distribution to lower column densities assuming the distribu¬ 
tion can be fit by a gamma function and constrained by the 
total number of LL systems. This is shown in Figure ^ by the 
dotted line. There is a tremendous drop (a factor of 100) in 
the number of absorbers after the cold gas truncation level. 
Of course a different source of absorbers may compensate 
for this drop, yet it seems unlikely that they would exactly 
cancel the drop caused by the truncation of cold gas disks. 
More data in this regime would be very helpful in clarifying 
the physical nature of the absorption systems. 

Since there does not seem to be enough cross section to 
produce the observed number of LL systems in either galac¬ 
tic halos or mini-halos, we must infer that the majority of 
LL systems are likely to be associated with uncollapsed gas 
in the vic inity of galactic h alos. This is consistent with the 


results of Dave et al 


(1999), who found that in a hydrody¬ 


namic CDM simulation, systems with column densities of 
~ 10^^ cm“^, typical of LL systems, are produced in regions 
with overdensities of ~ 100 at 2 ; = 3. One can see from these 
niilatinns that these systems tend to reside in the vreakly 
Qon-linear filaments surrounding halos. This gas may be en¬ 
riched by material ejected by supernovae from the halos and 
will probably have velocity widths ~ lOOkms”^. A line of 
sight passing through 1 kpc of 99% ionized gas with an av¬ 
erage over-density of 100 at 2 = 3 acquires enough optical 
depth to be above the Lyman limit. These are conditions 
that would not be unusual around galaxy mass halos. Fur¬ 
ther hydrodynamic simulations will be needed to study this 
Qon-linear regime in more detail. 

More data are required to identify the origin of LL ab¬ 
sorbers at high redshift and it is quite possible that hot gas 
in galactic halos, cold gas in mini-halos, and uncollapsed 
gas around halos all play a non-negligible role. Studies of 
the ionization state of the gas, its column density distribu¬ 
tion and its kinematics will help identify the nature of these 
absorption systems. 


5.3 The kinematics of sub-DLA systems 

A kinematic investigation of sub-DLA systems provides a 
generic, direct test of the multiple component model. In any 
multiple component model, most lines of sight pass though 
only a single component. The only way for DLA systems 
to be dominated by multiple components is for the column 
density of a single component to be below the DLA cutoff. 
This means that at some lower column density (sub-DLA) 
the absorption systems must become dominated by single 
component encounters. In our model this happens at HI 
column densities less than 10^°cm“^, as shown in Figure 
7. This figure shows the fraction of absorption systems pro¬ 
duced by more than one gas disk as a function of column 
density. The sub-DLA systems are very different below a log 
column density of 20. They are almost entirely composed of 
single disk systems. Thus we expect a significant change in 
the measured kinematics for these lower column density sys¬ 
tems. The difference in the measured An for the low ions is 
shown in the upper two panels of Figure ^ One sees the 
transition to many more small An systems. Observing the 
low-ion An in these systems should be a direct test of the 
multiple component model. 



log N 


HI 


log N, 


CIV 


Figure 7. 
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The fraction of systems composed of multiple gas disks 
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3ers. The DLA systems are mostly multiple compe¬ 
ls as required to explain their kinematics; however the 
(particularly with column densities below lO^'^cm”^) 
produced by single disks. CIV absorbers that are not 
with DLA or sub-DLA systems (dashed line) are more 
single component systems than CIV associated with 
ilid line). 
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xa et al. (2001) have studied a system at 2 = 2.5 
column density of 19.4. This system is neutral 
by our definition should be refered to as a sub- 
m even though its column density is less than the 
value we use in our model. The system has very 
ematics clearly indicating it is a single component 
iwever, it was also selected for its simple kinemat- 
y deuterium so it is unclear if it is representative 
'j ulation. An unbiased investigation of systems like 
needed to test the multiple-component scenario, 
s trong trend with column density seen in our model 
artifact of the fact that we have rather artificially 
all gas disks at the same value. A spread in val- 
tend to weaken the effect; nevertheless, the trend 
T for any multiple component model for DLA sys- 
use such a model must always have a larger cross¬ 
single encounters than multiple encounters. Thus 
the kinematics of sub-DLA systems is a critical 
multiple component model. 

liar effect would be predicted for CIV profiles if 
omposed of multiple components. Again, lines of 
gh single components are preferentially lower col- 
ty systems. The case of /sub = 1 is shown the right 
;ure 7 for CIV associated with DLA and sub DLA 
id for CIV with no associated low ions. The large 
in systems with column densities between 13.5 
uggests that it may be possible to identify which 
mtain cold neutral gas from the CIV kinematics. 
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Figure 8. The distribution of model velocity widths An is shown 
for the low ions of DLA systems, the low ions of sub-DLA systems 
and the high ions of DLA, sub-DLA and LL systems. While the 
high ions show little difference between DLA and sub-DLA sys¬ 
tems, the low ions are markedly different, providing a direct test 
of the multiple component model. A similar effect is seen between 
the high ions in the LL systems and the others, an additional test 
of our model. 


This provides yet another diagnostic of our model for the 
high ions. 

The difference between the velocity widths Avciv for 

f stems with and without low ions is quite striking (Figure 
lower three panels). This trend holds even if /sut, = 0 
where all systems are produced by a single component. In 
that case the DLA and sub-DLA systems arise from lines of 
sight passing through the central parts of the halo and thus 
are more likely to have large Avciv- Lines of sight that do 
not intersect cold gas are farther from the halo center and 
sample a smaller range of velocities. We note however that 
since some other state of gas is giving rise to most of the LL 
systems this gas may also produce CIV and thus complicate 
this picture. Further conclusions will not be possible without 
also modelling the other sources of LL absorption (i.e., the 
ones unaccounted for by our model). 


6 CONCLUSIONS 


We have presented a model for the high-ionization state gas 
seen in DLA systems, based on a galaxy formation model 
set within the hierarchical structure formation paradigm. 
The basis of this model is to associate hot halo or sub-halo 
gas with the highly ionized gas that gives rise to CIV ab¬ 
sorption. Assuming a simple model for the hot gas distribu¬ 
tion, we generate absorption systems with kinematic prop¬ 
erties in reasonable agreement with observations of DLAS 
at 2 = 3. We thus conclude that a CDM-based galaxy for¬ 
mation model can also account for both the low- and high- 
ionization state gas in observed DLA absorption systems. 
Agreement with the data is not a strong function of any free 
parameter of the model and thus we conclude that the gen¬ 
eral concept of a spherical hotter component as the source 
of CIV absorption provides a viable explanation for the ob¬ 
served kinematic properties. Associating some hot gas with 
subhalos and including a CIV component in co-rotation with 
the low ionization state gas improves the agreement, but 
models without these features cannot be categorically ruled 
out. 


While the galaxy formation model provides the infor¬ 
mation about the amount of hot and cold gas present in a 
given halo, our model requires a number of additional as¬ 
sumptions about the gas distribution. The most striking is 
our requirement that the cold gas be in a rather extended 
configuration (paper I). Also, in order for our model to work, 
we require that the gas termed “hot” in the semi-analytic 
mo del is in the form of c l ouds in a two-phase medium (as 
in Mo & Miralda-Escude 1996). This is clearly beyond the 
level of detail that semi-analytic galaxy formation models or 
hydrodynamic simulations currently attempt to address. 


Our model produces good agreement with the observed 
distribution of CIV column densities under the assumption 
that hot gas is distributed like the dark matter and that 
the fraction in a state suitable for CIV absorption fciv, 
and the metallicity, Zhg, are uniform throughout the halo. 
While this description has the advantage of simplicity, it 
is unlikely to be correct. The fact that the different high- 
ionization species do not always trace one another implies 
that fciv and Zhg vary throughout the gas. We found that 
a model in which a larger fraction of the gas was associated 
with sub-halos produced better agreement with the observa¬ 
tions, but an alternative scenario in which the subhalo gas 
possesses higher values of Zhg or fciv might provide a bet¬ 
ter explanation. The other high-ionization state metal lines 
may be used in future studies to help constrain the distribu¬ 
tions of gas density, metallicity and ionization state in high 
redshift halos. 


We also explored the predictions of our model for the 
lower column density sub-DLA and LL systems. We sug¬ 
gest that the kinematics of sub-DLA systems may be a de¬ 
cisive test of the multiple component model. Any multiple 
component model must have a large cross-section to single 
component intersections. Most sub-DLA systems should be 
produced by lines of sight passing through a single disk and 
therefore they should have significantly different kinemat¬ 
ics than the DLA systems. This will also be true for CIV 
systems if the hot gas is clumped into subhalos. Finally, we 
expect a large drop in the number of absorption systems 
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with column densities below the level where our gas disks 
are truncated, an additional test of the model. 

We find that in our model, only about a third of the 
observed population of LL systems are produced by hot gas 
in halos. We furthermore find in agreement with Abel & 
Mo ( 19|98| ) that there is enough cold gas in mini-halos to 
produce all of the LL systems at 2 = 3 and that there is also 
a similar amount of hot gas in these halos. However if we 
adopt a realistic model of the gas in mini-halos, we hnd that 
they produce an almost negligible contribution to the cross- 
section of LL systems. This is because to cover a large area 
the gas must have a low column density and therefore will 
be photoionzied by the UV background. Having high enough 
column densities to be self shielding implies the cross-section 
is small and therefore not a significant source of LL systems. 

We therefore propose that gas surrounding halos but 
outside the virial radius (e.g. in hlaments) may give rise to 
LL syste ms, a view support ed by the hydrodynamic simu¬ 
lations of Dave et al. (1999). If this gas is pre-enriched or 
enriched by metals ejected from the halo it may also pro¬ 
duce metal line systems. Further investigations using hydro- 
dynamical simulations will be useful in studying uncollapsed 
gas outside of virialized halos and determining its contribu¬ 
tion to absorption systems. 

We have illustrated that absorption systems comprise 
a powerful probe of galaxy formation. Acquiring more data 
of this kind, and developing more detailed models of these 
systems, will help in forming a complete picture of the prop¬ 
erties of the gas present in the early epochs of galaxy for¬ 
mation, which form the building blocks of the galaxies that 
we see today. 
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